Cold nuclear effects on heavy flavours (a review) 



R. Granier de Cassagnac 

8 ■ L8b " atoire & * P oly« e ch„i<,„ e /IN2P 3 . Palaiseau, 91128 France 

, E-mail: raphael@in2p3.fr 

(N ■ 

' Abstract. Before wondering about the quark-gluon plasma (QGP), one has to 

^) . take into account various cold (normal) nuclear matter effects, that can be probed 

through p+A like collisions. This article aims at reviewing the current results (and 
understanding) of these effects on heavy quarks and quarkonia production. 
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Qh| 1. Foreword 

The understanding of normal nuclear effects on heavy flavours production recently 
became crucial to understand the J/ip nuclear modification factors measured with 
^ ■ the PHENIX experiment [lj. In particular, the observed rapidity dependence of J/ip 
suppression is reversed with respect to what one would naively expect from density- 
t-h ■ induced suppression mechanisms, such as the colour screening originally proposed as a 
QGP signature [2]. J/ip are more suppressed at higher rapidity, where the density is 
lower. The question is: Could cold nuclear matter effects explain this behaviour? Or 
do we need to invoke the coalescence of c and c quarks coming from uncorrelated pairs 
instead? This question is so central that the two contributions on quarkonia suppression 
^ ■ in heavy ion collisions we heard at this conference discussed it in details [3j H]. In this 
paper, I focus on cold nuclear matter only, on most of the available p+A like collisions 
(from SPS, FNAL, HERA-B and RHIC) and on all the heavy flavour observables (J/i/j, 
ip', T and open charm...). In the following sections, I show the various dependencies of 
the heavy flavours productions, as they appeared chronologically. 

But before to look at data, one needs to define a somewhat arbitrary boundary 
between what we call "hot" and "cold". Initial state effects (EMC effect, parton 
antishadowing, shadowing or saturation, energy loss and Cronin effect, intrinsic 
charm...) obviously belong to the cold sector. Among the final state effects, absorption 
by incoming nucleons is usually classified as a cold effect too. On the contrary, 
absorption by comovers (of often unknown partonic or hadronic nature) is an effect 
of the hot matter. Obviously, so are the long awaited colour screening and quark 
recombination or coalescence. In a short and pragmatic way, I define as cold what 
can be grasped in p+A like collisions, while hot is what happens in addition when one 
looks at A+A collisions. In this paper, we shall then mostly look at p+A like collisions. 
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2. At SPS, some charmonia and a nuclear length scaling 

The left part of figure CD summarises the J/ ip over Drell-Yarfl productions in a large 
variety of collisions, as measured at the CERN SPS (y/SNN — 20 GeV) j7]. From 
p+p, to various p+A, to S+U, and up to peripheral In+In or Pb+Pb collisions, the 
data falls exponentially as a function of the nuclear thickness parameter L. The simple 
interpretation of this scaling is that J/ if) are absorbed by the forthcoming nucleons, the 
number of which being proportional to ij§|. From simultaneously taken p+A data, the 
NA50 experiment extract an absorption cross-section of 4.3 ± 0.5 mb for the J /if), as 
well as 7.7 ± 0.9 mb for the higher if>' state, as seen on the right part of figured] [8j. As 
expected, a lesser bound state is more fragile and suffers larger absorption. 

Beyond this normal absorption, J /if) are further suppressed in both the In+In [9] 
and Pb+Pb [TU] more central collisions (right-most points in figured] left), while the if)' 
already melts in S+U collisions [8]. 
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Figure 1. Left: J/ip yields normalized by Drell-Yan, as a function of the nuclear 
thickness L, as measured at the SPS. Right: J/tp and ip' yields in p+A collisions, 
normalised by A, and fitted to different absorption cross-sections. 



3. At FNAL and HERA-B, more quarkonia and xp scaling 

At slightly higher energy {Js^ ~ 40 GeV), the E866 (ELI and HERA-B p2] 
experiments have also measured that if)' were more absorbed than J /if). In the E866 case, 
this can be seen on the left part of figure [2] where the power a by which the production 

| As expected, Drell-Yan is not very affected by the medium, as it was reported by NA50 [5] and can 
also be seen on the right part of figure 2 from E722 [6|. 

§ Indeed, L is computed by counting the average number N of subsequent collisions in a Glauber model, 
converting it in a distance through L — N/crp, where p is the nucleon density and a the inelastic N+N 
cross section (which is finally arbitrary, cancelling out in the N / a ratio) . 
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is modified with respect to binary scaling (<j p a = <J PP x A a ) is depicted as a function 
of Feynman's refers to the momentum fraction of the parton in the 

proton projectile x p or target nucleus xa)- Focussing on the xf — region, we indeed 
see that ip' are more suppressed than J/ip. It is also the case in the negative Xp region 
(down to —0.35), probed by HERA-B [12]. In the higher xf region probed by E866, 
both the J I ip and ip' yields decrease dramatically, as shown on the left part of figure [2j 
This is not explainable without invoking effects beyond simple absorption cross-section, 
in particular gluon shadowing and energy loss of the incoming parton (see [13] or [2] 
for two complicated but successful fits of this data). However, since no A+A data exist 
in this region, we shall focus on xf — 0- 
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Figure 2. Quarkonia suppression measured at FNAL. Left: through the a parameter 
versus xf- Right: through the ratio of p+A to p+d cross-sections versus A. 

Other particles were significantly measured by these experiments. The upper 
diamond on the left part of figure [2] corresponds to D-mesons measured in E789 |15j . 
which (with an accuracy of ~ 4 %) do not seem to be modified by the medium (a ~ 1). 
This was also confirmed recently by HERA-B [IB]. The right part of figure [2] compares 
quasi-unmodified Drell-Yan [6] (diamonds) to suppressed charmonia [17] (circles) and 
bottomonia [18] (squares) in various p+A (normalised to p+d) collisions. Bottomonia 
suffer less than charmonia from the medium, indicating a smaller (effective) absorption 
cross-section. HERA-B recently released a detailed preprint [19] on Xc production, in 
which the ratio of J/ip coming from Xc lJ /ip is measured to be 18.8 ± 1.3 ± 2.4 %, 
not depending on the two (C and W) targets used for this measurement. 



4. A word on a possible a a b s energy dependence 

Having measured J/ip production at several energies, allows one to ask the question: 
does cr a bs vary with energy? While theoretical arguments exist to support the hypothesis 
that it decreases with energy, for instance because the Lorentz boosted and contracted 
nucleons are already gone when quarkonia form, no firm experimental conclusion exists 
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yet. The interplay between a a b s and (at least) shadowing, as well as the different x 
ranges probed by experiments make it difficult to perform such comparisons, as shown 
for instance in [20J. 

5. At RHIC, J/ip and their rapidity dependence 

At RHIC, the PHENIX experiment is able to measure J/ip at mid rapidity (\y\ < 0.35), 
but also at forward rapidity (1.2 < \y\ < 2.2), probing a large rapidity range. The 
left part of figure [3] shows the J/ip nuclear modification factor measured in d+Au 
collisions [2T] . Positive rapidity J/ip are originating from lower x partons in the Au 
nuclei (down to 2.10~ 3 in this case). They suffer more suppression, and this is usually 
interpreted as a sign of initial parton shadowing. At this energy (^snn = 200 GeV), 
cc pairs mostly originate from the fusion of gluons, the shadowing of which is not 
precisely constrained by existing data. Assuming various shadowing schemes, one can 
then derive the additional effect of nuclear absorption. As an example, on the left part 
of figure [3] assuming the NDSG [22] shadowing, a a a bs (called (?breakup here) value is 
extracted and quoted on the figure. The experiment confessed at this conference [3] 
that part of the systematic uncertainties was not properly propagated to a a bs- However, 
one can already see that cold nuclear effects at RHIC are not well constrained, the 
uncertainty on a a b s being at least of 2 mb. 
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Figure 3. J/ip nuclear modification factor in d+Au collisions as a function of rapidity, 
as measured by the PHENIX experiment. Left: Extraction of a break-up cross-section, 
assuming the NDSG scheme. Right: Illustration of the impact of intrinsic vs. extrinsic 
production mechanisms. 

More complications even arise from the underlying mechanism of J/ip production. 
Creating J/ip quasi-alone (inheriting its pt from the initial partons intrinsic p?) or with 
a hard gluon (balancing the J / ip extrinsic pt) should a priori probe different x values of 
the initial partons. Such an effect is derived from the EKS shadowing scheme [23J and 
compared to PHENIX data in a recent preprint |24j . Comparing the red open and closed 
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symbols on figure [3] (right) , which differ only by their intrinsic and extrinsic nature, the 
authors show that such an effect could change RdAu by about 10 %. This comes on top 
of a large uncertainty due to the shadowing scheme used. 

While it seems that nuclear absorption and shadowing are needed to explain J/ip 
production in d+Au collisions at RHIC, it is also clear that there is an interplay between 
the two, as it is for instance illustrated in [20], in which the authors derive a 1.7 mb 
systematic uncertainty on o a b s by varying the assumed shadowing scheme. It is also 
clear that gluon shadowing is poorly constrained by data, as the authors of the EKS 
model stress in their updated analysis leading to the new EPS scheme [25] . 

Two approaches are proposed to deal with this interplay and uncertainties. Trying 
to avoid relying on (shadowing) models will be addressed in the next section through 
a data driven method, while the status of open charm measurements, that would allow 
one to disentangle initial from final state effects will be shown afterwards. 

6. At RHIC, J/ip and their centrality dependence 

Unlike for lower energy p+A like collisions, d+Au centrality is measured at RHIC, and 
J /if) have been split in four centrality bins, at the three PHENIX rapidities. In |26j, 
it was proposed to take advantage of this impact parameter dependence to avoid 
relying on shadowing scheme. To do so, the author first performs phenomeno logical 
fits to RdAu(y,b), where b is the impact parameter derived from a Glauber model, 
also used to extrapolate to Au+Au. For each Au+Au collision occurring at a 
given impact parameter b^uAu, the N co u elementary nucleon-nucleon collisions are 
randomly distributed (following Woods-Saxon nuclear densities) providing the locations 
b\ and b\ of each collision i, relative to the centre of nucleus 1 and nucleus 2. For 
the considered Au+Au collision, the predicted nuclear modification factor RauAu is 
given by the following summation over the elementary collisions: RauAu(v, &a«a«) = 
X/iJ?" RdAu(—y, b\) x RdAu(+y, b % 2)/N co u. This method has the advantage not to depend 
on shadowing scheme, a a b s or production mechanisms, and to allow an easy propagation 
of the (large) RdAu experimental uncertainties to RauAu- The latest results of this 
calculation (from [21]) are shown on the left part of figure H] and compared to measured 
RauAu- The obtained authorised cold suppression (depicted as areas) have large enough 
uncertainties to allow the anomalous suppression to be equal at mid and forward 
rapidity. 

It is to be noted that this conclusion is also reached through the first attempt to 
derive J/ip production in the framework of the Colour Glass Condensate (CGC), in a 
recent preprint [27]. While the authors fit the absolute amount of J/ip suppression to 
the data itself, its rapidity dependence does not depend on this global fit. The right part 
of figure H] shows a comparison of their result with the measured RauAu- The bottom 
panel is the ratio RauAu(u = 1-7)/ RauAu(u = 0). While this CGC approach manages to 
reproduce the difference in central collisions, it fails in the more peripheral ones. 
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Figure 4. J ftp RauAu for the two rapidity measured by PHENIX, as a function of 
Npart, and compared to, left: an extrapolation of cold effects from the d+Au measured 
centrality dependence and to, right: a Colour Glass Condensate calculation. 

7. At RHIC, a poorly known open charm 

An interesting experimental way that could help disentangling J/ijj nuclear absorption 
and gluon shadowing would be to measure open charm with good precision. Open 
charm is indeed sharing its initial state effects (shadowing) with charmonia, while the 
final state effects should be totally different. At the moment, open charm is known to 
scale with the number of binary collisions, as we have seen on the left part of figure [2] for 
FNAL energies, and is shown in [28] for RHIC energies. The latest suffer a systematic 
uncertainty of the order of 25 % that can only be reduced with the help of silicon 
tracking detectors allowing to measure displaced vertices. 
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Figure 5. Nuclear modification factor of electrons from heavy flavour decays. Left: 
from STAR. Right: computed from PHENIX preliminary results. 

While the bulk of heavy flavour production scales with binary collisions, both 
PHENIX [28] and STAR [29J have reported a large suppression of electrons coming 
from heavy flavour decays at high p?- The situation in d+Au collision is less clear, as 
it can be seen on figure [5j The left part is the STAR measurement of the Au+Au (red 
circles) and d+Au (green squares) nuclear modification factors [25] , while the right part 
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reflects an extraction of the same quantity I made, based on preliminary d+Au [30J 
and published p+p [31J PHENIX data. Both figures exhibit an interesting enhancement 
that will be much better addressed by the analysis of the 2008 RHIC data which have 
30 times the statistics used here. If confirmed, such an enhancement could be interpreted 
in terms of antishadowing and/or Cronin effect, depending on its pt behaviour. 



8. And the Cronin effect everywhere 

The so-called Cronin effect consists in multiple scattering of an initial parton on the 
facing nucleus, resulting in an increased pt of the final state. In E866 for instance, a 
raise of a with pt is observed for Drell-Yan, charmonia [H] and bottomonia [18] . The 
average resulting p\ should vary proportionally to the amount of centres the partons can 
scatter upon, that can be characterised by L, the thickness parametejjj]. The left part of 
figure [6] shows < p\ > as a function of L, as measured by SPS experiments^!- For each 
energy, it clearly exhibits the expected linear dependence. The left part of figure [6] shows 
the current status of < p\ > measurements at RHIC [32J . I plot them here as a function 
of L and perform simple linear fits to all the points, including p+p, d+Au, Cu+Cu and 
Au+Au data. Within this limited accuracy, the mid rapidity slope is compatible with 
zero while the forward slope has a significance of 2.7 standard deviations, which could 
be a first sign of Cronin effect. This measurement clearly deserves more precision which 
should come soon for d+Au (2008 run) and Au+Au (2007 run). 
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Figure 6. J/tp <p\> versus L from (left) SPS experiments (right) PHENIX. 



9. LHC, the uncharted territory 

As a conclusion, let us think about the future LHC energy regime. In the previous 
sections, we saw that cold effects affecting heavy flavours were not well understood 

|| Indeed, for instantaneous processes, the average number of subsequent or preceding interactions are 
proportional, and the same L thickness parameter can be considered for initial (Cronin) of final (nuclear 
absorption) effects. 

% The missing NA60 p+A measurement was shown in conferences but not published. 
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when one considers a large amount of data (large xp coverage or various V /Sivjv). They 
are also not well constrained by experimental data, especially at RHIC. It would thus 
be very daring to predict them at LHC. As an exemple, J/vp at midrapidity will be 
produced from the fusion of two gluons with x < 10 -3 . The difference between various 
shadowing models in this regime is about a factor of two (to be squared in the fusion), 
for Q 2 ~ m 2 c . It will thus be crucial to have p+A like collisions at LHC, in order to 
interpret heavy flavour productions... 
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